Damping of long- wavelength collective excitations in quasi-onedimensional Fermi 

liquids 



F. Capurro, M. Polini, M. P. Tosi* 
INFM and Classe di Scienze, Scuola Normale Superiore, 1-56126 Pisa, Italy 



Abstract 

The imaginary part of the exchange-correlation kernel in the longitudinal current-current response function of a 
quasi-onedimensional Fermi liquid is evaluated by an approximate decoupling in the equation of motion for the current 
density, which accounts for processes of excitation of two particle-hole pairs. The two-pair spectrum determines the 
intrinsic damping rate of long-wavelength collective density fluctuations, which is calculated and contrasted with a 
result previously obtained for a clean Luttinger liquid. 
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1. INTRODUCTION 



An established theoretical viewpoint on one-dimensional (ID) fermion fluids is that the fermion-fermion interactions 
induce a Luttinger-liquid state, in which the momentum distribution does not have a discontinuity across the Fermi 
surface and spin-charge separation occurs [1]. A quasi-lD system which is of special interest both for basic physics 
and for technology is obtained by confining electrons in GaAs-based semiconductor quantum wire structures. The 
excitation spectrum of such a system in the ID quantum limit, where only the lowest subband is populated, has been 
measured by Goni et al. [2] by resonant inelastic light scattering. Focusing here on the intrasubband transitions, 
they consist of charge-density and spin-density collective excitations as well as of single-particle (electron-hole pair) 
excitations. As discussed by Wang et al. [3], no definitive evidence for Luttinger-liquid versus Fermi-liquid behaviour 
has been obtained from these data. The Fermi surface can be restored in a weakly interacting ID system by scattering 
of the electrons against impurities or thermal fluctuations [4] and at any rate the differences between a Luttinger liquid 
and a Fermi liquid are most marked in the one-electron spectrum, rather than in the particle-hole spectrum studied by 
Gofii et al. [2]. In particular, the plasmon dispersion relation that these authors have measured is well described by the 
standard random-phase (RPA) expression [5] , in a regime of wave number and wire thickness where the logarithmic 
factor from long-range Coulomb interactions is slowly varying so that an almost linear dispersion is observed. 

A number of authors have recently drawn attention to the possible relevance of the Tomonaga-Luttinger model 
in determining the density profile and the low- lying excitations of quasi- ID quantum gases of atoms inside highly 
anisotropic traps [6-8]. The extremely high purity and low temperature of these systems permit refined experimental 
studies using the tools of atomic physics and quantum optics. With regard to excitations in a clean Luttinger liquid, 
their lifetime has been calculated by Samokhin [9] in the case of short-range (delta-function) interactions. After 
stressing that a linear dispersion relation, for which the laws of energy and momentum conservation are simultaneously 
satisfied, leads in second-order perturbation theory to a divergent decay rate in the limit of an infinitely long sample, 
Samokhin carried out a self- consistent calculation of the damping of excitations from their mutual collisions that are 
allowed by a finite band curvature^ [10]. 

In the above perspective it seems useful to carry out a similar calculation for the lifetime of collective excitations 
in a ID Fermi liquid and to contrast the result with that obtained by Samokhin [9] for a ID Luttinger liquid. As is 
known from the study of interacting electron fluids in higher dimensionalities [11 - 13], the damping of a collective 
excitation at long wavelengths is determined by its decay into two particle-hole pairs and the rate of this process 
can be evaluated from the imaginary part of the exchange-correlation kernel entering the equation of motion for 
the current density. The same approach applied to a neutral Bose fluid describes the Beliaev-Popov mechanism for 
the damping of phonons [14] and more generally allows the evaluation of its visco-elastic spectra [15]. The only 
important conceptual difference for the quasi-lD liquid is that transverse currents are quenched in this case by the 
tight transverse confinement. For a quasi-lD Fermi liquid this damping mechanism has only been studied in the 
asymptotic high-frequency limit [16]. 

The paper is organized as follows. Section 2 explicitly shows how the inclusion of band curvature in a non- 
perturbative way leads to the same equations of motion for the particle density and current density in a quasi-lD 
Fermi fluid as for Fermi fluids in higher dimensionalities. We proceed from there to evaluate the exchange-correlation 
kernel at T = in Section 3 and the related damping rate of collective excitations in Section 4. Section 5 summarizes 
our main conclusions and comments on the role of temperature. 

2. THE ROLE OF BAND CURVATURE IN PARTICLE DENSITY FLUCTUATIONS 

We consider a clean ID liquid of interacting particles on a segment of length L ^ oo, which are described by the 
Hamiltonian [10] 
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Here and in the following spin indices and sums over them are left implicit, since wc are only interested in the total 
particle current density. In Eq. (1) p = ±1 is the branch index, while the operators c^.^ {ck,p) create (annihilate) a 
particle of momentum k and branch index p and obey the canonical fermionic anticommutation relations {cj, ^ , c^/ } = 
Sp,p'5k,k'- The symmetric density operator p{k) is defined as 



P Q,P 

and the single-particle spectrum Sa{pk) is taken as 

Sa{pk) = vpipk - kp) + {pk - kp)'^ , (3) 

with kp the Fermi momentum and vf = kp/m {h = 1). The linear Luttinger model (LLM) takes a = in Eq. (3), 
while the original Haldane spectrum is recovered by taking a = 1. 

We evaluate the collective particle-density excitations coming from the Hamiltonian (1) by the equation of motion 
approach, which was first introduced for the LLM by Metzner and Di Castro [17]. The relevant commutators are 

[pp{k), K] = Y^ [saipq) - £a{p{q - k))] ^-fe.pC^.p (4) 

Q 

and 

[p^{k),u]=pk^ m ■ (5) 

We separately treat below the two cases a = and a = 1. 



2.1 The LLM case (a = 0) 

For a = Eq. (4) yields [pp{k),K] = pvFkpp{k) and hence idtp{k) = vpkpas{k), where the asymmetric density 
operator is given by Pas{k) = pPCg_i. ^Cg^p. Using Eqs. (4) and (5) to get the equation of motion for pg^{k), one 
easily finds an undamped harmonic-oscillator equation for the density fluctuations, dfp{k, t) + Vi^{k)p{k, t) =Q where 

= (ypkf [1 + 2V{k)l{'KVF)] . (6) 

This is the result that is commonly obtained by the bosonization technique for the Luttinger liquid (see e.g. Samokhin 
[9]). It has also been obtained for a ID Coulomb fluid by the standard Bohm-Pincs RPA [18] and, taking in this 
case V{k) = 2e^Ko{kd) with Ko{x) the modified Bessel function and d the transverse width of the confinement, one 
recovers the well-known result fl{k) {^e'^VF/T^Y''^ k\ ln(fcd)|^/^ for the plasma frequency at long wavelengths. 

Starting from the undamped excitations described by Eq. (6), Samokhin [9] has evaluated the spontaneous decay 
of one excitation into two at temperature T = 0, as induced by the quadratic term in Eq. (3) which plays the role 
of an effective interaction between the Tomonaga-Luttinger bosons. His calculation assumes a constant value for the 
interaction potential (C/q, say) and yields from Eq. (6) excitations having a sound-wave dispersion relation, fl{k) = ck 
with c = VF[i- + 2Uo/{'!TVf)]^^'^- The imaginary part of the self-energy of the excitations at frequency ck has the form 
7fc oc at T = 0. 
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2.2 Current density in the case a = 1 



In the case a = 1 all linear terms cancel out inside the brackets in Eq. (4). One finds 

k 



[pp{k),K]=k ^_ 

and hence one obtains the continuity equation in the usual form, 



(7) 



idtp{k,t) = kJ{k,t) . (8) 
Here, the current density J{k, t) has the standard expression for a fluid of Fermi particles, 

Q,P Q,P 

(see e.g. Singwi and Tosi [19]). 

The equation of motion for the current density follows by taking the commutator of J{k) with the Hamiltonian (1), 

idJik, t) = —E{k, + 7 E -Vi<l)Pik - q, t)p{Q, t) (10) 

Q 

where the operator E{k) has been defined as 

^(^) = Yl ^4-fc/2,P^9+fc/2,P • (11) 
Q,P 



The equation of motion for particle density fluctuations follows from Eqs. (8) and (10), 



1 jL„ 

d^p{k,t) = -— E{k,t) - 7 E - V{q)p{k-q,t)p{q,t) . (12) 



The RHS of Eq. (12) contains both the restoring forces on density fluctuations and their damping by spontaneous 
decay and mutual scattering. The leading long-wavelength term for a ID Coulomb fluid is obtained in the RPA by 
preserving only the term gr = A; in the summation over momenta and by dropping the kinetic term: one flnds undamped 
oscillations at a frequency [nk'^V{k)/m]^/^ with n = N/L = 2A;ir/7r being the electron density. The limiting value is 
again the plasma frequency (4e^i;i?/7r)^/^A:| In {kd)}^^"^. 

In summary, we have seen that full account of the quadratic term in the single-particle spectrum in Eq. (3) formally 
leads to the same equations of motion for particle number and longitudinal current density fluctuations as for Fermi 
liquids in higher dimensionalities. We proceed in the next section to transcend the RPA by calculating the exchange- 
correlation spectrum for current density fluctuations at long wavelength. 



3. THE EXCHANGE-CORRELATION SPECTRUM 



The exchange-correlation kernel fxc\i^) for longitudinal currents in a fermion fluid is a quantity of crucial im- 
portance for time-dependent density functional theory [20, 21]. It is deflned as the long- wavelength limit of the 
function 



/W(A;,a;) = ^ 



1 



1 



V{k) , 



(13) 



x1{k,Lu)+n/m XL{k,uj) + n/m 

where XL(fc)'^) is the current-current linear response function and x"(^j'^) is the same function in the absence of 
interactions. ExpHcitly, the current response function is given by XL{k,uj) = (( J(fc); J(— fc)))^,, according to the 
general definition of a response function in terms of unequal-time commutators: 
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POO 

{{A; B))^ = -i dt {[A{ty, B{0)]) exp [z(a; + is)t] . (14) 
Jo 

Here A{t) is the operator in the Heisenberg representation, the brackets denote the ground-state expectation value, 
and £ is a positive infinitesimaL 

The last term on the RHS of Eq. (13) ensures that collective excitations do not contribute to the kernel at long 
wavelengths. In fact, to leading order in this limit the exchange-correlation spectrum is given by 
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WicH'^) = lim ^^ImxL(fc, w) (15) 

where XL{k,oj) is the proper current response function. The real part of the kernel can then be obtained from the 
spectrum by means of the Kramcrs-Kronig relation [11]. 

Returning to the expression for XL{k,uj) and following the treatment given by Nifosi et al. [11], we find that the 
current response function satisfies the equation of motion 

uj\L{k, u) = ([[J(fc), H\J{-k)]) - (([J(/e), H\- [J{-k), H]))^ (16) 

and hence a lengthy but straightforward calculation yields 

ImxL(A:, w) = V r(g, fc)r(g', -k)lm{{J{q)p{-q); J{q')p{-q')))^ + o{k^) (17) 

where 

r(g, k) = q\V{\q + fc|) - V{q)\ - qkV{q). (18) 

The four-point response function entering the RHS of the exact result in Eq. (17) describes the spectrum of particle- 
hole multipair excitations in the interacting Fermi liquid, with a leading contribution coming from both direct and 
exchange two-pair processes. 

Following again Nifosief al. [11], the direct contribution to the two-pair channel can be extracted from Eq. (17) by 
the approximate decoupling 

Im{{AB;CD))^ - - T — [Im((i; C')),Jm((B; + 

Jo 

+lm{{A;D))^,lm{{B;C))^-^,] . (19) 

Such a decoupling approximation should be corrected for exchange terms, which are known to be very small at low 
frequency both on the basis of a general argument [11] and of detailed calculations in higher dimensionalities [13]. 

However, the exchange corrections reduce the value of Im fxc\'-^) by a factor 1/2 in the high-frequency limit, which 
was studied in various dimensionalities by Glick and Long [22] and by Holas and Singwi [23]. 

By using Eq. (19) in Eq. (17) we obtain from Eq. (15) the exchange-correlation spectrum at long-wavelengths in 
the two-pair approximation in the form 

Im/(^)(u;) = -(27r2n2)-i g,{u;) / duj' / dq \V{q)\'' Imx(g,a;') Imx(g,a; - io'). (20) 

Jo J-oo 

Here, ai, is a numerical coefficient to be obtained from the value of the kernel in the high-frequency limit, and gx{to) 
is a factor accounting for pair exchange processes which interpolates between 1 at low oj and 0.5 at high w on their 
energy scale (twice the Fermi energy). In deriving the expression in Eq. (20) we have also made use of the exact 
relation 

XL(fc,a;)+n/m = a;\(fc,a;)/fc^ (21) 
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holding between XL{k,(-o) and the density response function x(fc,Cij). 

At the level of second-order perturbation theory, the fimctions x(fc,Cij) in the integrand in Eq. (20) should be 
replaced by their equivalents for the non-interacting Fermi fluid {xoikji^), say). In this case the integral diverges 
when V{k) is taken as the bare Coulonib interaction (the same divergence is in fact present for electron fluids in 
higher dimensionalities [11 - 13]). To account for all higher-order terms in an approximate way, one may draw on the 
analogy of the theory with the standard treatment of transport phenomena via a collision integral and Fermi's golden 
rule. This suggests that in the integrand in Eq. (20) one may replace x{k,Lv) by Xo{k,co) provided that at the same 
time the bare potential is replaced by an effective scattering potential. The use of Thomas-Fermi screening, or even 
the replacement of V{k) by a constant, are the simplest approximations that one may invoke for actual numerical 
calculations. 

In summary, the expression that we propose for the exchange-correlation kernel in the two-pair approximation reads 
Im/W(^) = _(27rV)-Ui 5,(a;) / dco' dq\Vsc{q)\''lmxo{q,uj')lmxo{q,uJ - uj') . (22) 

Jo J-oo 

Here, Vsc{k) is a screened Coulomb potential, that in the next section we shall take as a constant Uq to make connection 
with the experiments of Gofii et al. [2] on the one hand and with the damping of excitations in a Luttinger liquid 
calculated by Samokhin [9] on the other. 



4. LIFETIME OF EXCITATIONS 



From the work of Holas and Singwi [23] we find that in a ID Fermi liquid the high-frequency behaviour of the 
spectrum of the proper density response function from two-pair processes is given by 



m^kpk'^ 2 



(23) 



27r2(mtj)9/2 ° 

in the case of constant effective interactions. Equations (15) and (21) then yield 

The same asymptotic result is obtained from Eq. (22) provided gx{u)) -^1/2 and ul = 1/2. We have made use of the 
expression 



(24) 



Imxo(fc,a;) 



k mui 
2kp kpk 



1 1 



k rmv 
2kp kpk 



(25) 



for the spectrum of the ID ideal Fermi gas, with i?(x) being the Heaviside step function. 

We are now ready to evaluate the damping of collective excitations in the ID Fermi liquid, on the assumption that 
their dispersion relation is of the acoustic type. We write Wk = ck — ijk and for |7/;| <C ck we have [11] 



7. = -^Im/i^Hcfc). 



Equation (25) is again used to evaluate the integral in Eq. (22) at low frequency, with the result 

7fe = 777TK- Uok^. 



64k% 



(26) 



(27) 



In obtaining this expression we have taken (w) = 1 at low frequency. 
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5. CONCLUDING REMARKS 



In summary, we have calculated the decay of collective density excitations in a quasi- ID Fermi liquid from the 
process of decay into pairs of particle-hole excitations. We have assumed a short-range (contact) effective interaction 
and a linear dispersion relation. We have argued, however, that our final result in Eq. (27) is also approximately 
applicable to a quasi-lD Fermi liquid of electrons in semiconductor quantum wires, in a range of wave number and 
wire thickness where an approximately linear dispersion relation is observed for the plasmon in inelastic scattering 
experiments. 

We have found that within this model the dependence of the decay rate on wave number is of the type 7fc oc at 

T = 0. The temperature dependence of the lifetime can be calculated by means of the fluctuation-dissipation theorem 
[14] and one easily finds 7fe cx Tfc^ in the regime ksT 3> hck. Therefore, the behaviours of the decay rate within a 
Fermi-liquid model are clearly distinct from those predicted for a Luttinger liquid. For the latter Samokhin [9] has 
reported 7fe oc A;^ at T = and 7fe(T) a T^/^A;^/^ at high temperature. 
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